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As noticed in Ref. 1,, the Ore-Powell's classical calculation of the o-Ps^ 37 decay amplitude does not fulfill 
Low's theorem requirements for the low energy end of the photon spectrum. We reanalyze the implications of 
Low's theorem applied to this decay considering the interplay between the soft photon limit and the energy scales 
present in the e + e~ system. For energetic photons, the spectrum agrees with the Ore-Powell result, but deviates 
from it when the photon energy is comparable to the positronium binding energy. In this region it is found 
that bound states effects are essential to reach agreement with Low's theorem and can be accounted for in the 
framework of non-relativistic QED. 



1. Introduction 

The spin-triplet (orthopositronium: o-Ps) 37 
annihilation rate was first obtained by Ore and 
Powell |2] in the late forties and can be found 
in many textbooks. Although the leptons inside 
the positronium atom are bound by the Coulomb 
potential, the effects of such binding have been 
commonly neglected in the calculation of the o- 
Ps decay rate, being the process dominated by the 
short distance part where the leptons annihilate 
into three photons. The Coulomb force will be 
expected to affect the annihilation process appre- 
ciably only when a photon of small momentum is 
involved, as already quoted in the original paper 
of Ore and Powell. The analysis of the low energy 
photon spectrum of the oPs — > 37 decay made in 
Ref. ,3] takes into account these novel features by 
the use of the non-relativistic effective field the- 
ory methods developed to study bound states in 
QED and QCD gEEE]. 

The later work was motivated by the obser- 
vation of Pestieau and Smith pQ, that the Ore- 
Powell calculation of the o-Ps decay is in apparent 
contradiction with Low's theorem. While Low's 
theorem applied to this decay requires that the 
spectrum vanishes as E"^, with E-y the energy 
of the radiated photon, the standard calculation 
predicts a 0{E 1 ) behavior for E~ — > 0. We show 
that the Ore-Powell calculation can be reconciled 
with Low's theorem if the latter is carefully ap- 



plied considering all the energy scales present in 
the photon decay spectrum: the electron mass m, 
the binding energy of order ma 2 , and the hyper- 
fine splitting between the singlet and triplet states 
of order ma r. The Ore-Powell computation is 
valid for photon energies E 1 ^> ma 2 . When E 1 
is of order ma 2 , the Ps binding energy can not 
be neglected. The decay amplitude depends on a 
sum over an infinite set of excited Ps states, and 
can be written in terms of a Ps structure function 
which was computed in Ref. |3] using NRQED to 
leading order in the velocity (v) expansion. When 
E-y is of order ma 4 , the decay amplitude is dom- 
inated by the o-Ps — > p-Ps transition and it can 
be shown that the photon spectrum crosses over 
from an E-y behavior above the structure func- 
tion region to an E^ behavior below the p-Ps res- 
onance region. 

2. Low's theorem and o-Ps — > 37 decay 

Low's theorem jS] gives the amplitude for soft 
photon emission in the scattering of charged par- 
ticles. It states that the first two terms of the se- 
ries expansion in powers of the photon energy of 
a radiative amplitude X — > Yj may be obtained 
from a knowledge of the corresponding nonradia- 
tive amplitude X — » Y: 

e fi M» = ^+M 1 + G(k). (1) 
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One finds that A4q and Aii are independent of 
k and completely determined from the nonradia- 
tive amplitude To, its derivatives in physically al- 
lowed kinematic directions, and the electromag- 
netic properties of the particles involved .8, . 

The term A4q arises from the emission of a 
photon by ingoing or outgoing charged particles 
and is proportional to To times the universal fac- 
tor —Qi€ ■ Pi/k ■ pi, summed for all the external 
lines in the diagram. Note that if the nonradia- 
tive process involves no moving charged particles 
or is forbidden due to some selection rule, then 
Aio is identically zero. The term Aii can be ex- 
pressed as a function of the magnetic moments, 
the amplitude Tq and its derivatives with respect 
to internal variables which are not subject to con- 
straint, e.g. the energy and angles. 

Combining the amplitude behavior with that 
of the phase-space, the low-frequency form of the 
photon spectrum is 
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where A is proportional to |A4o| and B is the 
MqJ\A\ interference term. If Aio vanishes, the 
soft photon decay spectrum is of order T 7 dT 7 ; if 
both M. and yVfivanish, it is of order T 7 dT 7 . 

In the three-photon decay of o-Ps, one of the 
photons can have an arbitrarily small energy. The 
process can be then viewed as the radiative ver- 
sion of the o-Ps^ 27 decay. As the two-photon 
decay of o-Ps is not allowed by charge conjuga- 
tion invariance, the direct application of Low's 
theorem yields ./Wo,i = so that the o-Ps— > 37 
amplitude is of order 0(T 7 ), and the decay spec- 
trum is 



dr or 
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Et 
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as T 7 — > 0. This is in contradiction with the 



Ore-Powell spectrum, 
dr 37 2ma 6 " 



d.r; 



9tt 



x = ^ (4) 



which vanishes linearly with T 7 , as pointed out 
in Ref. [Tj. 

To understand the origin of the contradiction, 
it is worth noting that in the derivation of Low's 



theorem, one takes the limit T 7 — > and ne- 
glects all states other than those degenerate with 
the incoming and outgoing states, i.e. one uses 
T 7 <§; AT, where AT is the energy gap to ex- 
cited states. The amplitudes Aio,\ depend on 
the charge and magnetic moment couplings be- 
tween all the intermediate states degenerate with 
the initial or final states. A more general version 
of Low's theorem gives the decay spectrum for 
small T 7 without taking the strict T 7 — » limit. 
One treats all states with AT <C T 7 as degener- 
ate states, and includes them in the computation 
of charge and magnetic moment matrix elements 
for the purposes of Low's result. Which states are 
included in Low's theorem then depends on the 
magnitude of T 7 . 

In the case of Ps decay, consider the case where 
the photon energy is much larger than the binding 
energy. Then all Ps states (including o-Ps, p-Ps, 
radial excitations, etc.) are degenerate for the 
purposes of Low's theorem. In this case, there is a 
non-zero magnetic dipole matrix element between 
o-Ps and p-Ps, so that Ai\ does not vanish in 
this extended space of states. As a result, the 
decay spectrum vanishes linearly with T 7 . This 
is the approximation under which the Ore-Powell 
calculation is valid. For energies much smaller 
than the o-Ps-p-Ps hyperfine splitting, p-Ps as 
well as radial excitations are treated as excited 
states, the matrix element Ai\ vanishes, and the 
spectrum is of order T 7 . 

The discussion above could be also applied to 
the 37 annihilation of excited Ps states, like the 
1 Pi -> 37 decay 0. 

3. Orthopositronium Decay Amplitude 

The Ore-Powell annihilation amplitude can be 
obtained from the free-particle decay amplitude 
given, to lowest order in a, by the graph in 
Fig. ^ If the o-Ps momentum space wavefunc- 
tion is 4> (p) an< i the free e + e~ annihilation am- 
plitude is A(p), the bound-state decay amplitude 
to lowest order in v is 
d 3 p 



(27T>- 



-A( p )Mp) - A(0)M0) , 



(5) 



where V'o(O) is the o-Ps position space wavefunc- 
tion at the origin. Spin averaging \A\ and in- 
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Figure 1. Three photon annihilation graph. The 
graph is summed over the 3! permutations of the 
photons. 

tegrating over the three-body phase-space yields 
the differential decay rate written in Eq. (0J. 

The non-relativistic effective theory computa- 
tion (see details in Ref. 0) provides a system- 
atic way of including bound state effects in the 
o-Ps decay amplitude. The NRQED Hamilto- 
nian is constructed to reproduce the QED am- 
plitude when binding effects are neglected. Once 
the NRQED Hamiltonian has been determined, 
it can be used to compute the decay including 
binding corrections. 

The bound state dynamics is described by the 
Coulomb Hamiltonian for an e + e~ system in 
interaction with the quantized electromagnetic 
field: 
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Hq + H int 
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with x and p the relative position and relative 
momentum of the pair of leptons, and cr^, cr x 
the Pauli matrices acting on the electron and 
positron spinors. The electromagnetic interac- 
tions are through multipole interactions with the 
electric and magnetic fields and only the elec- 
tric and magnetic dipole interactions are shown 
in Eq. the higher multipoles being of higher 
order in v. 

The Coulomb Hamiltonian Hq is the leading 
term in the velocity power counting. The kinetic 
energy and Coulomb potential are of the same or- 
der in v. The energies and wavefunctions of Hq 
are thus the Coulomb wavefunctions with reduced 
mass m/2. The electric and magnetic dipole in- 
teraction terms are treated as perturbations. The 
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Figure 2. (a) Magnetic dipole graph for o-Ps an- 
nihilation. The solid square denotes the p-Ps an- 
nihilation vertex, (b) Electric dipole transition for 
o-Ps decay. The solid square denotes the n 3 P ,2 
annihilation vertex 

o-Ps decay amplitude has both short distance and 
long-distance contributions. The long distance 
part of the o-Ps decay amplitude entails the radi- 
ation of the soft photon (E 1 <C m) from the o-Ps 
bound state, whereas the short distance part is 



given by the on-shell e 



27 QED decay am- 



plitude where the photons are hard (E ~ m). 

The effective theory graphs which describe the 
o-Ps— ► 37 decay are shown in Fig. [21 The mag- 
netic term in H mt can induce a l 3 5i — > 1 1 Sq 
transition between Coulomb e + e~ states. The 
only allowed intermediate state in the dipole ap- 
proximation is the p-Ps ground state, with energy 
E p , so that E a — E p = AE^g. The electric dipole 
term — ex • E in H mi can change orbital angular 
momentum by one unit, allowing for transitions 
from ground state o-Ps to n 3 Po,2 states (n/ 1). 

To compute the o-Ps decay amplitude match- 
ing coefficient, one compares the effective theory 
result with the QED computation which neglects 
bound state effects. For the matching computa- 
tion, the energy separation between o-Ps and the 
intermediate states, E a — E n ~ 0(ma 2 ) is taken 
to be much smaller than the photon energy E-y. 

4. Results for the low energy spectrum 

The o-Ps— ► 37 decay spectrum from our effec- 
tive theory is written in terms of the magnetic 
and electric amplitudes obtained from the graphs 
in Fig. H[31, 
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ma 



\ 2 + Ik\ 2 



dx 9tt 
For the magnetic term a m , it is found 
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A-Ehfs = — ma 4 , T n = —ma 5 . 
12 ' p 4 

The electric amplitude is written in terms of the 
p-wave Coulomb Green's function py, which con- 
tains the sum over n 3 Po,2 states in Fig. 03: 



d|/i/ 4 Gi(0,»;K r )^ o (v).(10) 



7/ ^o(o) y 

A closed analytical formula for this integral has 
been given in Ref. ^U]. The magnetic and electric 
amplitudes behave as 



a m (Ery) 



a e (E 7 ) 



1 
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E 1 > AE hts , 
E 1 < AE 

hfs- 



2B, 



-Ely 3> ma 2 , 



i£ 7 <C ma 2 



(11) 



In the limit E 1 ^> ma 2 the Ore-Powell spectrum 
(Eq. (@J) is recovered, and it is shown that the 
sum of the magnetic and electric dipole transi- 
tions in the effective theory gives the full theory 
amplitude at leading order in the non-relativistic 
expansion. The matching condition, which is the 
difference of the two results, vanishes 3 and there 
is no additional three-photon annihilation term in 
the NRQED Hamiltonian. 

The ratio of the photon spectrum to the Ore- 
Powell value is shown in Fig.[3]up to E 7 ~ ma/2. 
At energies large compared with the binding en- 
ergy, a e and a m approach their values in Eq. Ijlljl , 
and the spectrum reproduces the QED one. The 
magnetic and electric dipole terms contribute in 
the ratio 3:7. At energies small compared with 
the binding energy, the electric dipole transitions 
decouple, and one is left with the magnetic term 
due to the p-Ps resonance contribution, of 3/10 of 
the Ore-Powell value. At energies much smaller 
than the hyperfme splitting, the p-Ps state also 
decouples, and the decay rate vanishes as E^, as 
predicted by Low's theorem. 

The results are consistent with the Ore-Powell 
spectrum and with Low's theorem. They include 
binding effects in a systematic expansion in pow- 
ers of v. 

Acknowledgements 

We wish to thank S. Narison and his team for 
the organization of the QCD 04 conference. This 



E (keV) 
0.001 0.01 0.1 




0.8 



0.6 



0.4 



0.2 - 



100 



E (Ryd) 



Figure 3. Ratio of the o-Ps decay spectrum in- 
cluding binding energy corrections to the Ore- 
Powell spectrum (1 Ryd= ma 2 /2). 

work has been supported in part by EU Con- 
tract No. HPRN-CT-2002-00311 (EURIDICE), 
by MCYT (Spain) under grant FPA2001-3031, 
and by EU ERDF funds. 

REFERENCES 

1. J. Pestieau and C. Smith, Phys. Lett. B524 
(2002) 395. 

2. A. Ore and J.L. Powell, Phys. Rev. 75 (1949) 
1696. 

3. A. V. Manohar and P. Ruiz-Femenia, Phys. 
Rev. D 69 (2004) 053003. 

4. W. E. Caswell and G. P. Lepage, Phys. Lett. 
B 167, 437 (1986). 

5. G. T. Bodwin, E. Braaten and G. P. Lepage, 
Phys. Rev. D 51, 1125 (1995) [Erratum-ibid. 
D 55, 5853 (1997)]. 

6. M. E. Luke, A. V. Manohar and I. Z. Roth- 
stein, Phys. Rev. D 61, 074025 (2000). 

7. A. Pineda and J. Soto, Nucl. Phys. Proc. 
Suppl. 64, 428 (1998). 

8. F. E. Low, Phys. Rev. 110, 974 (1958). 

9. D. Gromes, Eur. Phys. J. C 36 (2004) 169. 

10. M. B. Voloshin, Mod. Phys. Lett. A 19 (2004) 
181. 



4 



